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We compute the axial vector current, form factor, and coupling for quarks in the instanton liquid 
model with two light flavors. Non-local current corrections are derived, as required by the effective 
't Hooft interaction. We obtain a pion-mediated axial form factor and an axial coupling which, 
pq ' when simply applied in the non-relativistic limit for constituent quarks, matches the experimental 

, value to within a few percent, both in and out of the chiral limit. 

o 

p . 1. INTRODUCTION 

<D ' 

The axial coupling constant of the nucleon, deduced from measurements of neutron beta decay, is known to be 
\Ga/Gv\ = 1.2670 ± 0.0035 in vacuum. Strict chiral invariance would have this ratio be unity, and the substantial 
difference suggests a deviation in excess of current quark mass effects. Thus it has long been argued that the increase 
in Ga is a result of the spontaneous chiral symmetry breaking that characterizes the QCD vacuum. 

Although additional chiral-invariant meson-nucleon couplings can remedy this problem in the linear sigma model 
Q , a more fundamental explanation is warranted. The issue is complicated by the difficulties experienced by lattice 
QCD practitioners ||, in that such calculations tend to yield values significantly lower than experiment. Possible 
explanations have recently centered on finite-size effects Q] , although the specific physics is unresolved || . 

One might also try to deduce the nucleon's axial coupling by studying that of a constituent quark, denoted in this 
£NJ ■ paper as gA- In the non-relativistic limit, 

O ; 

^ ■ G A = -g A , (1) 

Oh: 3 

but incomplete knowledge of the quark wave function leaves the relativistic corrections unknown. However there has 
been considerable success in treating the nucleon as massive constituent quarks j(3| , and if the deviation of G a from 

yS^, • unity is in fact due to chiral forces the constituent approach should encompass the relevant physics. At the quark 
level chiral symmetry also demands an axial coupling of unity, which generates an excess of nucleon axial charge when 
Eq. (Q) is directly applied. Studies addressing this predict negative corrections of 0(1) to gA, using the MIT bag 

^ " model ]7j , in the large- N c limit j|, ||] , or via perturbative corrections ]l(J . 

In this paper we consider corrections to gA from instantons. The Instanton Liquid Model (ILM), in which the QCD 



vacuum is populated by classical gauge configurations 11 , 12J , dynamically generates the constituent quark mass and 



has consistently led to phcnomenological success when applied to problems related to spontaneous chiral symmetry 
breaking jl3| Q. As detailed below, we find a 25% reduction in the axial coupling from instantons, with only slight 
sensitivity to the chiral limit. This is in excellent agreement with experiment if corrections to Eq. (Q) are small. 

To determine the axial coupling we will first carefully analyze the axial current. Since instanton-induced interactions 
between quarks are non-local, the currents associated with conserved quantities are non-trivial. Next, given the 
necessity of pions in axial phenomena, meson correlation functions are briefly reviewed. We then derive the axial 
current and verify the importance of the pion pole in the axial form factor, h,A(q 2 )- Finally, the leading corrections to 
the coupling gA are computed, taking into account multi-instanton effects and substantial non-local vertex corrections. 
In the chiral limit we find gA = 0.73 and estimate that with a physical pion mass gA — 0.77. With naive application 
of Eq. (0), these correspond to a nucleon axial couping of 1.22 and 1.28, respectively. 



II. NON-LOCALITY AND CONSERVED CURRENTS 

The effects of an instanton liquid background are most economically encoded in an effective quark action [ p5| . 
Taking the exact zero-mode solutions to the Dirac equation and isolating them as the dominant low-energy effects 
leads to a 2A r /-quark vertex function jl6|], similar to the Nambu-Jona-Lasinio model of QCD jl?], [l8|, l1|. A well- 
studied simplification of this model, chiral Random Matrix Theory, has shown that chiral symmetry breaking via 
random overlap integrals reproduces Dirac eigenvalue correlations as computed on the lattice |^oj . Thus such an 



2 



approach seems to indeed contain the essence of non-perturbative QCD. Furthermore, the ILM relies on only two 
parameters: the instanton density and average instanton size. Both of these were phcnomenologically fixed long ago 
with the vacuum energy density and chiral condensate [|l5| |2l| and measured on the lattice (22|. These parameters 
determine the diluteness of the instanton liquid, the ratio of the average size to inter-instanton spacing, to be about 
1/3, a somewhat small parameter which allows for perturbative treatment of the instanton-induced vertex. 
The Eulcidean effective action for Nf flavors in the chiral limit is written [fL4[ , 



S 



- [ d 4 x ij>1 (x)i$ip(x) + i\ [ dU d 4 z Y[ [ ~ d 4 x f d 4 y f iJj}(x f )i^(x f - z, U)$(y f - z, U)iftip f (y f ) , (2) 

f 



where U is the instanton's 2 x N c color/spin orientation matrix and z is its position. The $(x) is the zero mode 
solution for fermions in the field of one instanton, and its Fourier transform is the form factor, 



f(p) = 2x 



h(x)K Q {x) - Io^K^x) + -h{x)Ki{x) 

x 



(3) 



x—pp/2 



which provides a natural momentum cutoff of scale p~ x . We have left out additional (quarkless) A terms which relate 
this coupling constant, which is in fact a Lagrange multiplier, to the instanton density, N/V. 

Because of non-locality Eq. (^|) is not invariant under symmetry transformations. Since bulk observables are often 
insensitive to this problem it is usually ignored, however when addressing currents the dependence is crucial. A 
literature exists in which non-local interactions arc modified to be invariant IBS 



29 and we follow 



^ £4 |2§ |2§ gg, |2S 

the same procedure here. Taking the non-local four-fermion interaction as the starting point, it then becomes a 
matter of multiplying each quark operator by a path-ordered exponential in the background of a source gauge field 



replacing the Euclidean quark fields as: 



%l>\x) 



W(z, x)ij){x), 
^{x)-f W{x,z)j . 



with the operator 



W(x, y) = Pexp 



-/ / ds^ap" 1 



As has been pointed out in the cited works, the operators W(x,y) must be path-ordered and the choice of path 
integrated over is not unique. However, we are primarily concerned with the longitudinal currents, for which results 
are independent of any particular choice of path |27j . 
Thus we write the modified action as 



S = - d 4 x V' 1 \x)i$ij>(x) 



N f 



+iX / dU d A z 



Y[ [- d 4 x f d 4 y f ipj(xf)"foW(x f ,z)j i^^(x f - z)&{y f - z)iftW(z,y f )ip j \y f ) 



f 

We now concentrate on the case of two light quarks, Nf = 2. 

The Noether currents can be evaluated directly from the action. In momentum space, 



(4) 



j a M 



ss 



We thus have, after Fourier transforming the fields and averaging over colors, 



-£p_,/,t 



(2n) 4 



(5) 
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where the left-handed non-local term is 

■ L ( ) _ A(2?rp) flf (giixi? _ JLgiigi* 

J^KQ) - n 2 -1 9192 I jl h N n jl 



,4 ,4 d 4 pi d 4 p 2 d 4 ki d 4 k 2 d 4 p' 
X {2n)* (2^) 4 (2tt) 4 (2tt) 4 (2tt) 4 

x r^e^-'/e-fa-^Vj^V+M-fe-fe) Ut( Pl ) 7 oT7ol 1^' x (*i)/(j/)/(fci) 

Jx V 1 J /i l i 

_ eix .( fel - P 0- i ,.( Pl + P2 -p'-fe)^. i(pi) [T ^( fel )]^ /(p!)/^')}^*^)^^)/^)/^) . (6) 

and the right-handed j^(q) is similar defined. Group indices have been suppressed; note that these are carried by the 
operators T. The instanton or anti-instanton at the core of each vertex naturally splits the quark fields into right- or 
left-handed spinors, ^Jr/l = |(1 ± 7s)^- Since we are only concerned here with spin structures of the identity matrix 
or 75, the decomposition is trivial. Carrying out the spatial integrations, we have 



q 2 N 2 -l 



5^5 



1 



f!9 2 I jl 32 ft 32 h 



d 4 pi d 4 p2 d 4 k\ d 4 k 2 
(2^)4 (2tt) 4 (2tt) 4 (2tt)* 



<KPl +P2 - ^1 



4(piboT7o , C' 1 (fa) [/(pi) - /(pi + <z)] /(fci) - 4/^ (pi) \ T Mki)T h /(pi) [/(fci) - f(h q)} 

X^ /2l2 (P2)^£ 2j2 (fc 2 )/(p 2 )/(fc 2 ). 



(7) 



In the event of spontaneous symmetry breaking, quark operators which overlap with the condensate can be paired 
into loops [ p6| |30f . When chiral symmetry is broken an effective quark mass is generated, dependent on four- momentum 

P = y/PWil as 

M(p) = Mf{pf . 

The constant M is the solution to a gap equation obtained by minimization with respect to the coupling A [|3l| , 



N 
V 



= 4N r 



d 4 p M{p) 2 
(2tt) 4 p 2 + M(p) 2 ' 



(8) 



and the constant N/V = (200 MeV) 4 is the instanton density. With the phenomenological p = 1/3 fm, one finds 
M = 345 MeV. In this case the current, Eq. M), reduces to 



3 L M 



+ 



q 

A(2ttp) 4 

Nr. 



d 4 P 
(27T) 



4 Vl(p) i M (P + ihoTjo - M{p)T 



f(p)f(p + q) 



d 4 k 
J2^Y 



f(k)- 



M(k) 



M{k) 



-Tr 



>^ L (p + q),{9) 



where the trace is over spin only. 

Further reduction of this still-unwieldy expression requires specifying a particular current and associated group 
elements T. The simplest example is the quark number current, B^, where T = 1: 



BM = J 



d 4 p 
(2^ 



As q — > 0, this becomes |M 



S M (0) 



^(P) \ 7m - »% [M(p + 9 ) - Af(p)] }> V^P + q) 



dM{p) 
dPn 



(10) 



^4(p) (in -i 



_d 4 p 
(2 



V'(p) • 



It is easy now to see that the (/-dependent term in Eq. (|lfj|), a result of the non-locality of the interaction, is necessary 
for current conservation. With Eq. ( ]To| ) and the Dirac equation, 



\jl>-iM{p)}^(p)=0, 



(11) 




+ 




+ 



FIG. 1: The series of instanton ladder diagrams summed to obtain meson correlation functions. 



we see the standard cancellation ensues and 



q fi B fi (q)=0. 



(12) 



This is the simplest illustration of the more general rule, which becomes more complicated when the axial isovector 
current is considered. But before analyzing this, we review the pion's derivation in the ILM. 



III. MESON CORRELATION FUNCTIONS 



The pion's role in axial charge currents has been long established. To be made explicit in the ILM, one must sum 
all ladder diagrams using the 't Hooft vertex, a series shown in Fig. [I]. This was originally done in Ref. |3^], although 
the advent of the effective action (Q) has since made this procedure simpler. 

In the chiral limit, one finds the expected pole in the pion propagator as q 2 — > 0. Furthermore, the use of non-local 
corrections to the axial currents has been shown to produce the required transversality in the mesonic current and 
properly relate the pion decay and renormalization constants Following the results of these cited works, we have 
the pion correlation function of Fig. fjj 



d 4 P 
5 — ( L 



{{^{p)r A l5^(p + q)) (^(p + q)T B ^{p))) 



The gap equation 



4iV 2 \iX N C M 2 
can be written |32| 

1 



d 4 p M (p)M (p + q)[p-(p + q)+ M (p)M (p + cff\ 
(2tt) 4 {p 2 + M{p) 2 ) \(p + q) 2 + M(p + q) 2 ] 



d 4 p M{pf 



iX N C M 2 J (2tt) 4 p 2 + M{pf 



(13) 



(14) 



and one finds, upon expanding in small q 2 



n 



AB 



^ ) W c q 2 I " r2 




d 4 p M{p) 2 ~ \pM{p)M'{p) + \p 2 M'(pf 
[p 2 + M{p) 2 f 



.9; 



11 xAB 



(15) 



The coupling constant is determined thusly by the compositeness condition |34j and, with the standard instanton 



parameters listed above, we find 



3.78. The pion and other channels were analyzed in an NJL model with 



alternative form factors by Plant and Birse, who found g^qq — 3.44 for the pion and the sigma parameters g aq q = 3.51 
and m a = 443 MeV ^7|. This required extracting the sigma pole, and instead of repeating this procedure or fitting the 
sigma correlation function ]35| we will simply scale our pion-quark coupling in a similar manner to obtain g G qq = 3.86. 
These two channels will be sufficient in the following analysis of the axial properties of the quark. 



IV. AXIAL CURRENT AND FORM FACTOR 



We now address the axial isovector current. A key experimental quantity is the matrix element of the current, 
^(.P2 — Pi), taken with initial and final state nucleons of momenta p\ and pi. Here we compute the analogous 
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quantity with quark fields: 

(V(pa)]A«(«jr)|^(pi)) = t»t(pa)y [ 7 „755A(g 2 ) + q^h A {q 2 )] u( Pl ) , (16) 

where the u'(p2) and it(pi) are the Euclidean spinor solutions of the Dirac equation for free quarks, and q = P2 — Pi- 
We are ultimately interested in the limit of pi = p 2 = 0, or on-shell quarks at zero momentum transfer, with the 
constant qa = <7a(0) the axial coupling. First we consider /^(g 2 ), the axial form factor. With a nonzero qa and 
current conservation, it is clear in the chiral limit that h,A(q 2 ) has a pole at q 2 — 0, identified as the pion. Through 
explicit construction, we will show this pole is manifest in the ILM. 

Inserting T — 75T/2 into Eq. (^), we have bare and single instanton (and anti-instanton) contributions in 

A W = / + / ^^( P )^ 5 ^M(p + q)+M(p) 

2,M2j;pf -Mf(j> + <l) I 7S4 /(*) [f (* + 1) + /(* - «)] m ^£L W + ^) • ( 17 ) 



iV c J ^ ^ • v J ( 27r )4 'wim 1 i ^ v */j k 2 + M ( fc2 ) 
Graphically, this is Fig. ||a. 

Clearly, no amount of labor will reveal this to be a conserved current; as the appearance of the chiral Goldstone 
boson is a multi-instanton effect, so too is axial vector current conservation. Pion tadpoles, as rendered in Figs. ^> 
and |^c, are necessary to remove the finite divergence of A" The first is a bare current insertion, the second the 
non-local piece. Their sum is 

(18) 

This cancels the final term in Eq. (|l7]) and the full current is 

A ^ q) = J ^^^{^ + ^^ M ^ + lM (p + ^}^Y 4 ' ip + q) - (19) 

The first term is the bare axial coupling, unity with the definition of Eq. ( |l6| ) , and the second is the pion pole carried 
by the axial form factor, 

h A (q 2 ) = -l M{P + q) 2 +M(P) . (20) 

Note that at long wavelength the leading divergent part of the form factor is the non-local, quark level version of the 
classic nucleon result, hA(q 2 ) — —2iMj^/q 2 . The pion pole similarly couples separately to each constituent quark. 
It is now clear, with the Dirac equation (11), that 



Equivalently, with the propagator, 



the Ward-Takahashi Identity is satisfied, 



S[p) 



q^(q)=0. (21) 
i> + iM[p) 



p 2 +M{pY 



Q^U = ~ Y b 5S( P + ^ + ^rSs] . (22) 




FIG. 2: Diagrams contributing to the axial form factor, /ia(? )• Circled crosses are current insertions and the dashed line 
denotes a pion. 

V. AXIAL COUPLING CONSTANT 

The axial form factor involved the effects of multiple instantons, with pions transferring i-channel momentum. 
Similar propagation in the s-channel contributes to the axial coupling constant. While still dominant, the pion is not 
the only relevant resonance in this channel. Since the sigma's coupling constant is of comparable size, and its mass is 
only slightly greater than that of a constituent quark, its effects will also be taken into account. Higher resonances, 
such as the p and A\ vectors, are both substantially heavier and have quark-meson couplings a mere third of the 
scalar and pseudoscalar [35|. With the leading contributions of order gf these are safely ignored. 

In Eq. (|l^) we have seen that the bare gA is unity. To this we first add the pion contribution of Fig. ||a, which is 

3a., _ 9l qq f d*p M(p) p 2 + 2M(p) 2 
9a 2M J (2tt) 4 p 1 [p 2 + M(p) 2 ] 2 ' 1 j 

After evaluation, the identical sigma contribution differs only by the coupling constant and massive propagator in the 
integrand, 

3a., _ 9l qq f d±p M(p) p 2 + 2M(p) 2 
9a ~ 2M J (2tt) 4 p 2 + m 2 a [p 2 + M(p) 2 } 2 ' 1 ' 

These expressions evaluate to —0.102 and —0.0425, respectively when we take the parameters described in Section 



[II 



The non-local contributions couple the axial current directly to the four-quark vertex, as shown in Fig. |3|c. With 
an internal pion, we find 



and for the sigma, 



3M_^ W f d 4 p M'(p) 
9a AM J (2tt) 4 p[p 2 + M(p) 2 } ' 1 ' 



3b.a = dqq f _d*P pM' (p) 

9a AM J (2n^(p 2 + ml)[p 2 + M(p) 2 }- 1 ' 



Numerically, we obtain —0.101 and —0.0236 for these corrections. Note that both numbers are comparable to those 
of Fig. ^|a, i.e. the non-local corrections are not suppressed in any systematic way. 
In the chiral limit, we therefore find 

9A = l~ 0-102 - 0.0425 - 0.101 - 0.0236 = 0.73 , 

in excellent agreement with the non-relativistically deduced value. 

The vector current, to which the axial coupling is compared experimentally, receives no contributions from the 
non-local vertices. Furthermore, a critical sign difference in the numerator as compared to Eqs. ( p3| ) and ( p4| ) leads 
to a finite but very small contribution from each diagram, summing to only two percent of the bare coupling. Thus 
we can ignore such instanton effects in the coupling gy and take it to be unity. 
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/ \ / \ 
^ (8) * <& ^ 

(a) (b) 



FIG. 3: Diagrams contributing to the axial coupling constant, §a- Circled crosses are current insertions and dashed lines 
denote a either a pion or sigma meson propagator. 

In order to estimate the correction from finite current quark masses, we assume larger masses for the internal meson 
propagators. With the physical pion mass, m n = 138 MeV, the instanton calculations of Hutter J35| generate a sigma 
mass of m a = 540 MeV. Reevaluating, we naturally find a deduction in these subtractions, 

g A = l- 0.079 - 0.036 - 0.094 - 0.021 = 0.77 . 

However, this is a mere 5% increase compared to the chiral limit, and in fact places the result extremely close to the 
pseudo-experimental value. 

We are ultimately interested in the nucleon axial coupling, and in the non-relativistic limit we find 

_ 5 _ f 1.22 for m n = , > 

UA ~ 3 9A ~ \ 1-28 for = 140 MeV ' 1 

The physical value, Ga = 1-2670 [Q 1 , lies in between. 

VI. CONCLUSIONS 

Through an analysis of current conservation in the instanton liquid model, we found that terms arising from the 
non-locality of the 't Hooft vertex were central to maintaining the symmetries of the theory. This can be thought 
of as a correction for the truncation of quark basis states, for we have isolated the zero modes and removed them 
from the dynamics. Reintroducing the symmetries to the volume of this extended interaction vertex via path-ordered 
exponentials recovers the symmetries, as has been noted by other authors some time ago |24| ^5|. A general form of 
current corrections has been computed here, followed by specialization to the axial iso-vector case. 

After these modifications the pion pole appears as a t-channel quark-antiquark correlation, governing each con- 
stituent quark's axial form factor and maintaining current conservation. Similar multi-instanton corrections, when 
applied with s-channel chiral mesons, lead to an axial coupling constant of the quark of gA — 0.73. In the non- 
relativistic limit this corresponds to a nucleon coupling of Ga — 1.22, within four percent of the experimental value. 
Away from the chiral limit, when the pion is set to its physical mass, agreement improves to within one percent. These 
numbers are obtained with the two standard inputs of the instanton vacuum, the number density and average size, 
and no additional parameters. While there are additional corrections from vector mesons, the significantly lower cou- 
pling constants and higher masses render them negligible. However, the non-local current terms lead to contributions 
comparable to the local ones and cannot be disregarded. 

Our results suggest that the ILM contains the primary effects which lead to the observed vacuum axial coupling of 
1.267. Given that this variance from unity is a likely result of spontaneous chiral symmetry breaking, a central result 
of the ILM, it is not surprising that instantons contribute favorably. It is, however, noteworthy in its implication that 
rclativistic wave function corrections to the constituent quark are nearly negligible. 



Our convention differs from that reported by the Particle Data Group by an overall sign. 
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